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(Dated: February 18, 2001) 
We present measurements of time-dependent CP-violating asymmetries in neutral B decays 
to several CP eigenstates. The measurement uses a data sample of 23 million T(4S') — > BB 
decays collected by the BABAR detector at the PEP-II asymmetric B Factory at SLAC. In this 
sample, we find events where one neutral B meson is fully reconstructed in a CP eigenstate 
containing charmonium and the flavor of the other neutral B meson is determined from its 
decay products. The amplitude of the CP-violating asymmetry, which in the Standard Model 
is proportional to sin2/3, is derived from the decay time distributions in such events. The result 
is sin2/3 = 0.34 ± 0.20 (stat) ± 0.05 (syst). 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



CP-violating asymmetries in the time distributions for 
decays of 5° and mesons provide a direct test of the 
Standard Model of electroweak interactions ||l|. For the 
neutral B decay modes reported here, corrections to CP- 
violating effects from strong interactions are absent, in 
contrast to the modes in which CP violation was 
discovered 

Using a data sample of 23 million BB pairs recorded 
at the T(4iS') resonance by the BABAR detector at the 
PEP-II asymmetric-energy e'^e~ collider at the Stanford 
Linear Accelerator Center, we have fully reconstructed 
a sample Bcp of neutral B mesons decaying to the CP 
eigenstates J/ijjK^, tp{2S)Kg and Jf-ipK^. We examine 
each of the events in this sample for evidence that the 
other neutral B meson decayed as a 5° or a desig- 
nated as a P° or P° flavor tag. The final Bcp sample 
contains about 360 signal events. 

When the T(4S') decays, the P-wave BB state evolves 
coherently until one of the mesons decays. In one of four 
time-order and flavor configurations, if the tagging meson 
Ptag decays first, and as a P*^, the other meson must be a 
B° at that same time itag ■ It then evolves independently, 
and can decay into a CP eigenstate Bcp at a later time 
tcp- The time between the two decays At = tcp — 
ttag is a signed quantity made measurable by producing 



the T(4S') with a boost /?7 = 0.56 along the collision 
(z) axis, with nominal energies of 9.0 and 3.1 GeV for 
the electron and positron beams. The measured distance 
Az w P^cAt between the two decay vertices provides a 
good estimate of the corresponding time interval At; the 
average value of |Az| is PjcTgo « 250 /im. 

The decay-time distribution for events with a P*' or a 
P° tag can be expressed in terms of a complex param- 
eter A that depends on both P^P*^ mixing and on the 
amplitudes describing P" and P'^ decay to a common fi- 
nal state / . The distribution f+ (f_ ) of the decay rate 
when the tagging meson is a B^{B^) is given by 



f±(At) 



HAtl/T^o 



2rBo(l-H|A|2) 



i + |Ar 



± Xm\ sin (Am^o At) ^ 



cos (Ato^o At) 



, (1) 



where Tga is the B^ lifetime and Am^o is the mass dif- 
ference determined from P°P° mixing Q , and where the 
lifetime difference between neutral P mass eigenstates is 
assumed to be negligible. The first oscillatory term in 
Eq. is due to interference between direct decay and de- 
cay after mixing. A difference between the B^ and P*^ 
distributions or a At asymmetry for either tag is evidence 
for CP violation. 
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If all amplitudes contributing to i?" — > / have the same 
weak phase, a condition satisfied in the Standard Model 
for charmonium-containing b ccs decays, then |A| = 1. 
For these CP eigenstates the Standard Model predicts 
A — rife~^'^^ , where r/f is the CP eigenvalue of the state / 
and /? = arg [ -^dKb/^td^tb ] is an angle of the Unitar- 
ity Triangle of the three-generation Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [^. Thus, the time-dependent 
CP-violating asymmetry is 



Acp{At) 



f+{At) - f_(At) 
f+(At)-|-f_(At) 
— ry/ sin2/3 sin (Am^o At), 



(2) 



where 77/ = -1 for J/ijK° and ■>jj{2S)K° and -1-1 for 

A measurement of Acp requires determination of the 
experimental At resolution and the fraction of events in 
which the tag assignment is incorrect. A mistag fraction 
w reduces the observed asymmetry by a factor (1 — 2w). 

Several samples of fully reconstructed mesons 
are used in this measurement. The Bcp sample con- 
tains candidates reconstructed in the CP eigenstates 
J/tPK° {K", -> 7r+7r",7r07rO), ^i2S)K° {K° ^ n+n-) 
and J/ipK^. The J/'i/' and ^/'(25') mesons are re- 
constructed through their decays to e^e~ and 
the ^{2S) is also reconstructed through its decay to 
J/?/'7r+7r~. A sample of B decays -Bfiav |§ used in the 
determination of the mistag fractions and At resolu- 
tion functions consists of the channels = 
7r+,p+,a^') and J/ijK*^ {K*° K+tt'). A control 
sample of charged B mesons decaying to the final states 
J/t/>X(*)+, 'ij;{2S)K+ and D'^*^°ii+ is used for validation 
studies. 

A description of the BABAR detector can be found in 
Ref. [0. Charged particles are detected and their mo- 
menta measured by a combination of a silicon vertex 
tracker (SVT) consisting of five double-sided layers and a 
central drift chamber (DCH), in a 1.5-T solenoidal field. 
The average vertex resolution in the z direction is 70 /ini 
for a fully reconstructed B meson. We identify Icptons 
and hadrons with measurements from all detector sys- 
tems, including the energy loss (dE/dx) in the DCH and 
SVT. Electrons and photons are identified by a Csl elec- 
tromagnetic calorimeter (EMC). Muons are identified in 
the instrumented flux return (IFR). A Cherenkov ring 
imaging detector (DIRC) covering the central region, to- 
gether with the (]E/dx information, provides K-tt sepa- 
ration of at least three standard deviations for B decay 
products with momentum greater than 250 MeV/c in the 
laboratory. 

We select events with a minimum of three recon- 
structed charged tracks, each having a laboratory polar 
angle between 0.41 and 2.54 rad and impact parameter 
in the plane transverse to the beam less than 1.5 cm from 
the beamline. The event must have a total measured en- 
ergy in the laboratory greater than 4.5 GeV within the 



fiducial regions for charged tracks and neutral clusters. 
To help reject continuum background, the second Fox- 
Wolfram moment 1^ must be less than 0.5. 

An electron candidate must have a ratio of calorimeter 
energy to track momentum, an EMC cluster shape, a 
DCH dE/da; and a DIRC Cherenkov angle (if available) 
consistent with an electron. 

A muon candidate must satisfy requirements on the 
measured and expected number of interaction lengths 
penetrated, the position match between the extrapolated 
DCH track and IFR hits, and the average and spread of 
the number of IFR hits per layer. 

A track is identified as a kaon candidate by means of 
a neural network that uses dB/dcc measurements in the 
DCH and SVT, and comparison of the observed pattern 
of detected photons in the DIRC with that expected for 
kaon and pion hypotheses. 

Candidates for J/ip £^£^ must have at least one 
decay product identified as a lepton (electron or muon) 
candidate or, if outside the calorimeter acceptance, must 
have DCH dE/dx information consistent with the elec- 
tron hypothesis. Tracks in which the electron has radi- 
ated are combined with bremsstrahlung photons, recon- 
structed as clusters with more than 30 MeV lying within 
35 mr in polar angle and 50 mr in azimuth of the pro- 
jected photon position on the EMC. The second track of 
a fJ-'^ IJ-" pair, if within the acceptance of the calorime- 
ter, must be consistent with being a minimum ionizing 
particle. Two identified electron or muon candidates are 
required for J/t/j or ip{2S) reconstruction in the 

higher-background '0(2S')ifs and J/tpKl channels. 

We require a J/4> candidate to have 2.95 < m^+^- < 
3.14GeV/c2 or 3.06 < m^,+^,- < 3.14GeV/c2, and a 
ij{2S) -» i'^l" candidate to have 3.44 < m^+e- < 
3.74GeV/c2 or 3.64 < m^+^- < 3.74GeV/c2. Require- 
ments are made on the lepton helicity angle in order to 
provide further discrimination against background. For 
the '>p{2S) J/-07r+7r~ mode, mass-constrained J/ip can- 
didates are combined with pairs of oppositely-charged 
tracks considered as pions; the resulting mass must be 
within 15MeV/c2 of the ^(25) mass ||. 

A K'^ — ^ TT+TT^ candidate must satisfy 489 < TOtt+tt- < 
507MeV/c^. The distance between the J/^ or 4'{2S) and 
K^: vertices is required to be at least 1 mm. 

Pairs of 7r° candidates with total energy above 800 MeV 
are considered as K^^ candidates for the J/ipK^ mode. 
We determine the most probable decay point along 
the path defined by the initial A'!? momentum vector and 
the J/tp vertex by maximizing the product of probabili- 
ties for the daughter tt*^ mass-constrained fits. Allowing 
for vertex resolution, we require the displacement from 
the J/ip vertex to the decay point to be between —10 and 
-1-40 cm and the tt^tt^ mass evaluated at this point to be 
between 470 and 550MeV/c^ 

A candidate is formed from a cluster not matched 
to a reconstructed track. For the EMC the cluster must 
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FIG. 1: a) Distribution of jties and A_E for Bcp candi- 
dates having a Ks in the final state; b) distribution of AE 
for J/tjj K1 candidates. 



have energy above 200 MeV, while for the IFR the cluster 
must have at least two layers. We determine the en- 
ergy by combining its direction with the reconstructed 
J/ip momentum, assuming the decay B° J/ipK°. 
To reduce photon backgrounds, EMC clusters consistent 
with a tt" 77 decay are rejected and the transverse 
missing momentum of the event projected on the 
candidate direction must be consistent with the mo- 
mentum. In addition, the center-of-mass J/ip momentum 
is required to be greater than 1.4GeV/c. 

Bcp candidates used in the analysis are selected by 
requiring that the difference AE between the energy of 
the Be p candidate and the beam energy in the center-of- 
mass frame be less than three standard deviations from 
zero and that, for K^^ modes, the beam-energy substi- 
tuted mass mES = \/(-^bIfam)^ ~ (Pb")^ ff^lls in the inter- 
val rriES > 5.2GeV/c^. The resolution for AE is about 
10 MeV, except for J/i^Kl (3 MeV) and the 
mode (33 MeV). For the purpose of determining num- 
bers of events, purities and efficiencies, a signal region 
mES > 5.27GeV/c^ is used for all modes except J/-0i4r°. 

Figure H shows the resulting AE and tties distributions 
for Bcp candidates containing a K^, and AE for the can- 
didates containing a K^. The Bcp sample is composed of 
890 events in the signal region, with an estimated back- 
ground of 260 events, predominantly in the J/ip chan- 
nel. For that channel, the composition, effective rjf and 
AE distributions of the individual background sources 
are taken either from a Monte Carlo simulation (for B 



decays to J/tp) or the mg+g- sidebands in data. 

For flavor tagging, we exploit information from the in- 
completely reconstructed other B decay in the event. 
The charge of energetic electrons and muons from 
semileptonic B decays, kaons, soft pions from D* decays, 
and high momentum charged particles is correlated with 
the flavor of the decaying b quark: e.g., a positive lepton 
yields a B^ tag. Each event is assigned to one of four 
hierarchical, mutually exclusive tagging categories or is 
excluded from further analysis. The mistag fractions and 
efficiencies of all categories are determined from data. 

A lepton tag requires an electron or muon candidate 
with a center-of-mass momentum pcm > 1-0 or 1.1 GeV/c, 
respectively. This efficiently selects primary leptons and 
reduces contamination due to oppositely-charged leptons 
from semileptonic charm decays. Events meeting these 
criteria are assigned to the Lepton category unless the 
lepton charge and the net charge of all kaon candidates 
indicate opposite tags. Events without a lepton tag but 
with a non-zero net kaon charge are assigned to the Kaon 
category. 

All remaining events are passed to a neural network 
algorithm whose main inputs are the momentum and 
charge of the track with the maximum center-of-mass mo- 
mentum, and the outputs of secondary networks, trained 
with Monte Carlo samples to identify primary leptons, 
kaons, and soft pions. Based on the output of the neural 
network algorithm, events are tagged as B^ or 5° and as- 
signed to the NTl (more certain tags) or NT2 (less certain 
tags) category, or not tagged at all. The tagging power 
of the NTl and NT2 categories arises primarily from soft 
pions and from recovering unidentified isolated primary 
electrons and muons. 

Table | shows the number of tagged events and the sig- 
nal purity, determined from fits to the toes (^s modes) 
or AE {K1 mode) distributions. The measured efficien- 
cies for the four tagging categories are summarized in 
Table H 

The uncertainty in the At measurement is dominated 
by the measurement of the position ztag of the tagging 
vertex. The tagging vertex is determined by fitting the 
tracks not belonging to the Bcp (or Sflav) candidate to a 
common vertex. Reconstructed K'^ and A candidates are 
used as input to the fit in place of their daughters. Tracks 
from 7 conversions are excluded from the fit. To reduce 
contributions from charm decay, which bias the vertex 
estimation, the track with the largest vertex contribu- 
tion greater than 6 is removed and the fit is redone until 
no track fails the requirement or fewer than two tracks 
remain. The average resolution for Az — zcp — ztag is 
190 /im. The time interval At between the two B decays 
is then determined from the Az measurement, includ- 
ing an event-by-event correction for the direction of the 
B with respect to the z direction in the T(4S') frame. 
An accepted candidate must have a converged fit for the 
Bcp and -Stag vertices, an error of less than 400 /xm on 
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TABLE I; Number of tagged events, signal purity and result 
of fitting for CP asymmetries in the full CP sample and in 
various subsamples, as well as in the -Bflav and charged B 
control samples. Purity is the fitted number of signal events 
divided by the total number of events in the AE and tties 
signal region defined in the text. Errors are statistical only. 



Sample 




Purity (%) 


sin2/? 




273 


96 ± 1 


0.25 ± 0.22 




256 


39 ±6 


0.87 ± 0.51 


Full CP sample 


529 


69 ±2 


0.34 ± 0.20 


J/^K'^,^{2S)K"s only 


J/ipK% {K'i Tr+TT") 


188 


98 ± 1 


0.25 ± 0.26 




41 


85 ±6 


-0.05 ± 0.66 




) 44 


97 ±3 


0.40 ± 0.50 


Lepton tags 


34 


99 ±2 


0.07 ± 0.43 


Kaon tags 


156 


96 ±2 


0.40 ± 0.29 


NTl tags 


28 


97 ±3 


-0.03 ± 0.67 


NT2 tags 


55 


96 ±3 


0.09 ± 0.76 


B" tags 


141 


96 ±2 


0.24 ± 0.31 


W tags 


132 


97 ±2 


0.25 ± 0.30 


Bflav sample 


4637 


86 ± 1 


0.03 ± 0.05 


Charged B sample 


5165 


90 ± 1 


0.02 ± 0.05 



TABLE If: Average mistag fractions Wi and mistag diff'er- 
ences Awi = Wi{B^) — ■Wi{B^) extracted for each tagging 
category i from the maximum-likelihood fit to the time dis- 
tribution for the fully-reconstructed B" sample (Baav+Bcp)- 
The figure of merit for tagging is the effective tagging effi- 
ciency Qi — ei{l — 2wif , where Si is the fraction of events 
with a reconstructed tag vertex that are assigned to the i"* 
category. Uncertainties are statistical only. The statistical 
error on sin2/3 is proportional to l/\/Q, where Q ~^Qi- 



Category 


e{%) 


w (%) 


Aw (%) 


Q (%) 


Lepton 


10.9 ± 0.4 


11.6 ±2.0 


3.1 ± 3.1 


6.4 ± 0.7 


Kaon 


36.5 ± 0.7 


17.1 ± 1.3 


-1.9 ± 1.9 


15.8 ± 1.3 


NTl 


7.7 ± 0.4 


21.2 ±2.9 


7.8 ± 4.2 


2.6 ± 0.5 


NT2 


13.7 ± 0.5 


31.7 ± 2.6 


-4.7 ± 3.5 


1.8 ± 0.5 


All 


68.9 ± 1.0 






26.7 ± 1.6 



A2: and a measured | Az| < 3 mm; 86% of the Bcp events 
satisfy this requirement. 

The sin2/3 measurement is made with an unbinned 
maximum likelihood fit to the At distribution of the com- 
bined Bcp and iJflav tagged samples. The At distribu- 
tion of the former is given by Eq. |l|, with |A| = 1. The 
latter evolves according to the known rate for flavor os- 
cillations in neutral B mesons. The amplitudes for Bcp 
asymmetries and for -Bflav flavor oscillations are reduced 
by the same factor (1 — 2w) due to mistags. The distri- 
butions are both convoluted with a common At resolu- 
tion function, and are corrected for backgrounds, incor- 
porated with different assumptions about their At evolu- 
tion and convoluted with a separate resolution function. 
Events are assigned signal and background probabilities 
based on fits to tties (all modes except J/ipKl) or AE 



{J/ipKl) distributions. 

The At resolution function for signal candidates is rep- 
resented by a sum of three Gaussian distributions with 
different means and widths. For the core and tail Gaus- 
sians, the widths are scaled by the event-by-event mea- 
surement error derived from the vertex fits; the combined 
rms is l.lps. A separate offset for the core distribution 
is allowed for each tagging category to account for small 
shifts caused by inclusion of residual charm decay prod- 
ucts in the tag vertex; a common offset is used for the 
tail component. The third Gaussian (of fixed 8ps width) 
accounts for the fewer than 1% of events with incorrectly 
reconstructed vertices. Identical resolution function pa- 
rameters are used for all modes, since the -Stag vertex 
precision dominates the At resolution. 

A total of 35 parameters are varied in the final fit, in- 
cluding the values of sin2/3 (1), the average mistag frac- 
tion w and the difference Aw between B° and i3° mistags 
for each tagging category (8), parameters for the signal 
At resolution (9), and parameters for background time 
dependence (6), At resolution (3) and mistag fractions 
(8). The determination of the mistag fractions and signal 
At resolution function is dominated by the high-statistics 
Bflav sample, while background parameters are governed 
by events with toes < 5.27GeV/c^ (except J/i^K^). We 
fix tbo = 1.548ps and Amso = 0A72hps-'^ [|. The 
largest correlation between sm2(3 and any linear combi- 
nation of the other free parameters is 0.076. 

The measurement of sm2(3 was performed as a blind 
analysis by hiding the value of sin2/? obtained from the 
fit, as well as the CP asymmetry in the At distribution, 
until the analysis was complete. This allowed us to study 
statistical and systematic errors without knowing the nu- 
merical value of sin2/?. 

The measured mistag rates obtained from the likeli- 
hood fit for the four tagging categories are summarized 
in Table |l[ As a check, the mistag rates were evaluated 
with a sample of about 16,000 D*~ t^vi events and found 
to be consistent with the results from the hadronic decay 
sample. 

The combined fit to the CP decay modes and the flavor 
decay modes yields 

sin2/3 = 0.34 ± 0.20 (stat) ± 0.05 (syst). 

The decay asymmetry Acp as a function of At and the 
log likelihood as a function of sin2/3 are shown in Fig. ||. 
If |A| is allowed to float in the flt, the value obtained is 
consistent with 1 and there is no significant difference in 
the value of —rjfIm\/\X\ (identified with sin2/3 in the 
Standard Model) and our quoted result. Repeating the 
fit with all parameters fixed to their determined values 
except sin2/3, we find a total contribution of ±0.02 to the 
error on sin2/3 is due to the combined statistical uncer- 
tainties in mistag rates, At resolution and background 
parameters. 
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FIG. 2: The raw asymmetry in the number of and B^^ tags 
in the signal region, {Ngo — N-go)/{Ngo + Ngo), with asym- 
metric binomial errors, as a function of At for a) the J/tP Kg 
and 'tp{2S)K^ modes {rjf = -1) and b) the J/ip Kl mode 
{rif = +1). The solid curves represent the time-dependent 
asymmetries determined for the central values of sin2/3 from 
the fits for these samples. Eight events that lie outside the 
plotted interval were also used in the fits. The probability of 
obtaining a lower likelihood, evaluated using a Monte Carlo 
technique, is 60%. c) Variation of the log likelihood as a func- 
tion of sin2/3 for the modes containing Kg (dashed curve), the 
J/ip K1 mode (dotted) and the entire sample (solid) . For the 
latter, solid lines indicate the central value and values of the 
log likelihood corresponding to one statistical standard devi- 
ation. 



The dominant sources of systematic error are the as- 
sumed parameterization of the At resolution function 
(0.04), due in part to residual uncertainties in the SVT 
alignment, and uncertainties in the level, composition, 
and CP asymmetry of the background in the selected 
CP events (0.02). The systematic errors from uncertain- 
ties in Atobo and r^o and from the parameterization of 
the background in the selected i?flav sample are found to 
be negligible. An increase of 0.02 ?ips^^ in the assumed 



value for Am^o decreases sin2/3 by 0.012. 

The large sample of reconstructed events allows a num- 
ber of consistency checks, including separation of the 
data by decay mode, tagging category and i?tag flavor. 
The results of flts to these subsamples are shown in Ta- 
ble I for the high-purity events. Table | also shows 
results of fits with the samples of non-CP decay modes, 
where no statistically significant asymmetry is found. 

Our measurement of sin2/3 is consistent with, but im- 
proves substantially on the precision of, previous deter- 
minations The central value is consistent with the 
range implied by measurements and theoretical estimates 
of the magnitudes of CKM matrix elements jl0| ; it is also 
consistent with no CP asymmetry at the 1.1a level. 
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